Variational approximations to homoclinic snaking 
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We investigate the snaking of localised patterns, seen in numerous physical applications, using 
a variational approximation. This method naturally introduces the exponentially small terms re- 
sponsible for the snaking structure, that are not accessible via standard multiple-scales asymptotic 
techniques. We obtain the symmetric snaking solutions and the asymmetric 'ladder' states, and also 
predict the stability of the localised states. The resulting approximate formulas for the width of the 
snaking region show good agreement with numerical results. 
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There has been much recent interest in the phe- 
nomenon of spatially localised patterns, 0,0,0, Ell: ex- 
tending our understanding of earlier work on this topic 
[H, [13] ■ As discussed in the review article by Dawes , 
this work is motivated by wide-ranging applications in 
many different areas of physics, including buckling of 
struts and cylinders (1, [lol llHj , nonlinear optics Q , con- 
vection patterns [H [I4I , gas discharge systems and gran- 
ular media. Most theoretical work concentrates on the 
Swift-Hohenberg equation, which is the simplest model 
equation that illustrates the effect [l|, 0, |ll[ . 

These localised patterns arise as a result of bistability 
between a uniform state and regular patterned state. A 
front linking these two states might be expected to move 
in one direction or the other, except at a specific param- 
eter value, referred to as the Maxwell point. However, 
due to a pinning effect first described qualitatively by 
Pomeau [l3|, the front can lock to the pattern, resulting 
in a finite range of parameter values around the Maxwell 
point where a stationary front can exist. Combining two 
such fronts leads to a pinned localised state. The pin- 
ning range is seen in numerical simulations 14. 18| 
and leads to a 'snaking' bifurcation diagram in which the 
control parameter oscillates about the Maxwell point as 
the size of the localised pattern increases. 

The localised states can also been seen from the view- 
point of spatial dynamics as a homoclinic connection 
from the uniform state to itself, leading to a geometri- 
cal argument for the existence of such states over a finite 
range of parameter values d, H, [13] ■ 

It has been appreciated for some time [l^j that the pin- 
ning effect cannot be described by conventional multiple- 
scales asymptotics, since this method treats the scale of 
the pattern and the scale of its envelope as independent 
variables, leading to an arbitrary phase in the envelope 
function. Hence, the pinning range is 'beyond all orders', 
or is exponentially small in the small parameter corre- 
sponding to the pattern amplitude [l[. A very thorough 
analysis of the exponential asymptotics of this problem 
has recently been carried out b y K ozyreff and Chapman 
0, [HI , extending earlier work 1^ 3] ■ The calculation 
of the pinning range is extremely complicated and unfor- 
tunately requires two fitting parameters. 

Many previous authors have noted the variational 



property of the Swift-Hohenberg equation [2, [3, [l2| , but 
have not made use of this in regard to the snaking di- 
agram. Wadee and Bassom [3] did use the variational 
method, but not in the snaking regime. In this Letter we 
exploit the variational structure of the Swift-Hohenberg 
equation in order to approximate the pinning range. An 
ansatz based on the weakly nonlinear solution is substi- 
tuted into the Lagrangian of the system, and then this is 
minimised over the unknown parameters. This approach 
has a number of advantages: the calculations are much 
less cumbersome than those involved in the full exponen- 
tial asymptotics analysis; exponentially small terms ap- 
pear naturally through the evaluation of the Lagrangian 
integral; there are no unknown constants that have to be 
determined by fitting with numerical results; and the cor- 
responding effective Lagrangian can predict the stability 
of the states. Although the method is widely applica- 
ble, for illustrative purposes we consider the cubic-quintic 
Swift-Hohenberg equation [3, [3, [l3l • The results are then 
compared with numerics. 

The governing equation is given by 



dtu — ru — [1 + u + I 



(1) 



where r is the control parameter and and are the 
coefficients of the cubic and quintic term, respectively. 
The uniform state m = is unstable for r > 0. The 
interesting case leading to snaking bifurcations is 63 > 0, 
&5 > 0, in which case localised solutions can exist for 
r < 0. It is possible to rescale u to set 65 = 1 but we will 
keep 65 as a parameter for consistency with [2, [3] • 

The stationary solutions of ([Ij can be derived from a 
Lagrangian 
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which is finite since we are concerned with localised so- 
lutions. Furthermore, 

dC _ 



(dtuf dx 



(3) 



so stable solutions correspond to minima of C and C can 
be interpreted as the total 'energy' of a solution. The L^- 

norm of the solution is defined as iV = I v? dx\ 
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Motivated by weakly nonlinear analysis [H 3, we 
study two forms of localised pattern. For small r, we 
write the localised solution of ([1]) as 



— ^sech(i?x) cos(fca; + (p). 



(4) 



Using, e.g., complex variable techniques, substituting 
the approximation ^ into the Lagrangian ^ yields the 
effective Lagrangian 
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where 



5A%) +5A%k^ 
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Ki = 56B^ + ^A^B"^ (-81 

80B^ (fc2 - 1) + 45" {-10A%3 + 5A\ 

20A:^ + 30 (1 - r)) , 
K2 = AA^ {B^ + 4:k^) {B^ (-563 
K3 = A%5 {AB^ + AbB^k^ + 81/c'*) 
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Note that the terms involving the phase shift y in (O 
have an exponential factor. Therefore, it is expected that 
the splitting of localized solutions with different phase in 
the limit r — t' 0, i.e. — >■ 1 and i3 — >■ 0, is exponentially 
small, as suggested in [? , 11, 1^ . For the ansatz (|4]) 
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which is the square of the norm. 

Applying the Euler-Lagrange formulation to the effec- 
tive Lagrangian 



eff 



dsCeff = dkCeff = d^Ceff = 0, (6) 



gives us a system of nonlinear equations for A, B,k, and 
(fi that make (|4]) an approximate solution of ([TJ. Ne- 
glecting the exponentially small terms, ([6]) can be solved 
perturbatively to yield (up to O {{—r)^^'^)) 
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Note that the leading order expansions above are the 
same as those obtained using multiple scale expansions 
The phase-shift at this order is arbitrary, which 



also agrees with the multiple scales result. However, tak- 
ing into account the equation d^pCeff = 0, in which all 
the terms are exponentially small, ip is a multiple of n/2, 
i.e. there are two different types of solution, odd or even 
in X. This same conclusion was reported in after a 
lengthy 'beyond all orders' calculation. 

It is expected that the ansatz (jlj only resembles ex- 
act solutions of the Swift-Hohenberg equation when 63 
is order 1 and r is small. Computations of snaking are 
generally carried out with neither r nor bs being small 
[2,13, [3 ■ localised states in the snaking regime, the 
envelope solution has a plateau, which becomes longer as 
the norm N increases. As the Maxwell point is the center 
of the snaking, it is reasonable to construct homoclinic 
snaking solutions using a special solution at that point, 
i.e. a front solution. 

We approximate the solution in the snaking region by 



Acos(fca; + (p) 
Vl + e-B(l^|-^' 



(10) 



When the modulus sign in the denominator is absent, 
the solution ([TU| forms a front solution at the Maxwell 
point obtained using a multiple-scales expansion method 
for small b^ and small r [l|, [1, • Due to the modulus 
sign, ((TU)) patches two fronts of opposite polarity, lead- 
ing to a localised state of length 2L. Even though the 
derivative of this solution is not continuous at a; = 0, it 
is smooth enough when BL ^ 1, which we assume to be 
the case (i.e. the fronts are well separated). To simplify 
the algebra we set fc = 1. 

Substituting the ansatz pH)) into the Lagrangian 
the effective Lagrangian can be obtained as 
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+A^B^n{12 - B^ ~ 32r)] sin(2L) cos(2^) 
+ ^e-T [-30^A'^55 + 327rA% 
+ttB^A^{B^ - 4)] cos(2L) cos(2(p) 
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To illustrate some of the steps in the calculation, consider 

the term f °° v?dx. This can be written as 

J —00 
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t^dx — — — 
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after using the double-angle formula and the fact that 
the envelope function is even. The first term in the in- 
tegral can be evaluated directly and gives the answer 
A^L + 0{e~^^). The second term can be found by 
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writing cos(2x) =Re(exp(2ia;)) and using contour inte- 
gration. The integrand has simple poles at a: = ±L + 
in/B, ±L + Sin/B, . . ., but the first of these dominates, 
giving an exponentially small contribution of order e~"S" . 
Evaluating the residue in the standard way gives the 
leading-order result 



u'^dx = A^L 



2 

27r — er^ sin(2L) cos(20). (12) 
B 



The other terms in C can be evaluated in a similar way, 
using poles of order up to three. The formula (fT2|) gives 
the square of the norm N , indicating that in the snaking 
region, solutions with a larger norm correspond to longer 
plateaus. 

Consider first the terms in (fTTj) that are not exponen- 
tially small, and assume that _B ^ 1 so that the B^ term 
may be neglected. Making these terms stationary with 
respect lo L, A and B and solving them for A, B and r 
gives A = Am, a = aM, and r = tm, with 



Am = 3 



63 , 27bl 



(13) 



which are the parameter values of the front at the 
Maxwell point, in exact agreement with the multiple- 
scales asymptotic method P, 0, [l3l ■ Note that Cef / de- 
pends linearly on L, but this dependence vanishes at the 
Maxwell point. This is to be expected since the Maxwell 
point is determined by the condition that the patterned 
state has the same energy as the zero state. 
Near the Maxwell point we set 

r = tm + Sr. 

After making the above simplifications, 
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the patterned region grows in this case and shrinks if 
5r < 0. 

Finally we have the exponentially small terms respon- 
sible for the snaking. The first of these, proportional to 
sin2L, is larger than the second, for small 63. 

The two remaining undetermined parameters in the 
ansatz ([T0| are (f and L. The variation of the effective 
Lagrangian with respect to L and ip can be readily cal- 
culated from ([T4| . From d^Ceff — 0, there are two dif- 
ferent possibilities. The first is that sin(2(/3) = 0, so as 
in the previous case, (^c is a multiple of 7r/2. Solving the 
equation dhC^jf — for 5r at leading order in 63 then 
gives 



[448O65 + 576^] sin(2i) 
(14) 



In this form, the Lagrangian can be interpreted more 
easily. The first two terms are constant and independent 
of L. These arise from the two front regions at the ends 
of the localised states. The first of these terms dominates 
for small 63, being of order 6|. Since this term is positive, 
localised states have a greater energy than the periodic 
state or the zero state, so energy considerations suggest 
that the periodic state or the zero state are 'preferred' 
over localised states. Similarly, multi-pulse localised so- 
lutions will have a larger value of C, in proportion to the 
number of pulses. 

The third term in (1141) . proportional to —L5r, indicates 
a preference for increasing values of L ii Sr > 0, so that 
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(15) 



the plus and minus signs corresponding to the odd and 
even states respectively. This in turn gives the maximum 
value of Sr as 



Srr, 



IO65 



(16) 



The width of the snaking region is therefore approxi- 
mately given by 2Srm- 

The Lagrangian is then an oscillating function of L, 
intertwining between the Lagrangian with ip = and 
if = tt/2. There is a periodic sequence of alternating 
stable and unstable equilibria, corresponding to minima 
and maxima of the Lagrangian respectively. As Sr is in- 
creased, all the equilibria disappear in saddle-node bifur- 
cations at Sr — Srjn- Near these bifurcations, the stable 
state is the one of the pair with the lower value of L. 

The second possibility arising from d^Ceff = is that 
sin(2L) — (considering only the larger of the two ex- 
ponentially small terms in (fT4)) ). These types of solu- 
tions correspond to the 'bridge' 17 1 or 'ladder' [1] states 
that link the snaking branches, and only exist for val- 
ues of L that are multiples of 7r/2. For these solutions, 
dh^eff = determines the value of the phase by 

560&57raM . ,_,„^ 

Sr — ± —2 e cos(2(^). (17) 

903 

Considering Cef / as a function of the two variables (p 
and L, it is easy to see that the snaking solutions are 
either maxima or minima, while the ladder states are 
always saddle points. Therefore, the snaking solutions 
are either stable, or unstable with two positive eigenval- 
ues, and the ladder states are always unstable with one 
positive eigenvalue. These results are in agreement with 
numerical simulations (see [1]). 

We have found steady localized states of ([ij numer- 
ically, using a pseudo-arclength continuation method 
with periodic boundary conditions, implemented with a 
Fourier spectral discretisation. A summary of the results 
is shown in Fig. [T] 
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FIG. 1: (Top) The bifurcation diagram obtained numerically 
and our approximation obtained from solving ((6| for 63 = 
2. (Middle) The numerically computed Lagrangian ((2)1 as a 
function of the length of the plateau 2L, for 63 = 1. The 
dashed line is our approximation calculated from the effective 
Lagrangian (Bottom) The width of the snaking region 

as a function of 63. Filled circles are numerical and solid 
lines are analytical, i.e. 2Srm (|16p . The inset shows the same 
comparison in a log scale. In all the figures, — 1. 



The top panel shows the bifurcation diagram showing 
two branches of localized solutions for 63 = 2 and — 1. 
In the same panel, shown in dashed lines are our analyt- 
ical results ([6]), showing that the variational calculation 
approximates the numerics better for relatively small |r|. 

In the middle panel, we plot the Lagrangian ([2]) for 
&3 = 1 as a function of the length of the solution's plateau 



2L, which is calculated numerically as 2L w J v? dx, 

where um = max{u} and the integration is a definite 
integration over the computational domain. Plotted in 
the same panel is our effective Lagrangian (jlll) . There 
is good agreement for the average numerical value of 
the Lagrangian and the qualitative nature of the os- 
cillations, but the variational approximation underesti- 
mates the amplitude of the oscillations. The amplitude 
of the oscillations increases with L since from with 
6r cx cos(2L), there are oscillations of the form Lcos(2i). 

In the bottom panel we show the width of the snaking 
region as a function of 63 numerically and analytically, 
where our approximation is in fairly good agreement. 

To conclude, we have used variational methods to 
study the snaking behaviour of localised patterns in the 
Swift-Hohenberg equation. The approach has several ad- 
vantages: the exponentially small terms responsible for 
the phase locking arise naturally in the Lagrangian, giv- 
ing a simple formula p4l) from which the snake and lad- 
der localised states can easily be found, along with their 
stability. We have concentrated on a simple model equa- 
tion in the small-parameter regime, but the method is 
very widely applicable and can be expected to open up a 
new avenue of research into this challenging field. 
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